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Abstract A method is presented for finding anisotropic distribution functions for stellar systems with known, 
spherically symmetric, densities, which depends only on the two classical integrals of the energy and the mag- 
nitude of the angular momentum. It requires the density to be expressed as a sum of products of functions of 
the potential and of the radial coordinate. The solution corresponding to this type of density is in turn a sum of 
products of functions of the energy and of the magnitude of the angular momentum. The products of the density 
and its radial and transverse velocity dispersions can be also expressed as a sum of products of functions of the 
potential and of the radial coordinate. Several examples are given, including some of new anisotropic distribu- 
tion functions. This device can be extended further to the related problem of finding two-integral distribution 
functions for axisymmetric galaxies. 

Keywords celestial mechanics • stellar dynamics • galaxies 



1 Introduction 

There is a long history of studying the structure of galaxies by construction of self-consistent distribution func- 
tions for a stellar system with a known gravitational potential. The potential of the system determines the self- 
consistent mass density p of the system via Poisson's equation generated by the well-known Newtonian gravita- 
tional law, and also the structure of the stellar orbits according to Newton's equations of motion. The system is 
constructed from the building blocks of the orbits that can lie within the potential. This is called the "from p to 
/" approach for finding a self-consistent distribution function / (Binney and Tremaine 1987, hereafter BT). The 
distribution function (hereafter DF) of the system represents how stars are distributed in the phase space of the 
system and the integration of the DF over velocity space yields the density. Therefore the problem of finding the 
DF is that of solving an integral equation; in the spherical case in particular this equation is of the first kind, in 
which the unknown DF / only occurs inside the integral. 

Some outstanding astronomers have contributed to this mathematical problem. Eddington (1916) showed 
that it can be solved for an isotropic DF that depends only on the energy in the spherical case, whose velocity 
dispersions in the radial and tangential directions must be equal, by first expressing the density as a function 
only of the potential, and then solving the Abel integral equation. It is also well known that the isotropic DF is 
unique for the potential of any spherical system. 

Fricke's (1952) expansion method can be obviously generalized to the spherical case (e.g., Camm 1952; 
Bouvier 1962, 1963; Ossipkov 1979a; Kent and Gunn 1982; Dejonghe 1986, 1987; Dejonghe and Merritt 1988; 
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Cuddeford 1991; Louis 1993), and similarly yields the result that DFs which are products of the two powers of 
the energy and the square of the angular momentum correspond to densities which are proportional to products 
of the potential and the spherical radial coordinate for spherical systems. Hence the DF for the system can be 
obtained by first expressing the density as a function of the potential and the spherical radial coordinate, and 
then expanding as a power series (Qian and Hunter 1995). Moreover, there may be an infinity of anisotropic DFs 
corresponding to any given mass density in spherical stellar systems (Dejonghe 1987). 

A class of typical anisotropic DFs that depend on the energy and the magnitude of the angular momentum, 
whose velocity dispersions are anisotropic (but ellipsoidal), was independently found by Ossipkov (1979) and 
Merritt (1985), for a spherical density distribution. Such anisotropic DFs of Ossipkov-Merritt type are in fact 
analogues of Eddington's formula, mentioned above. In a number of papers, different integral transformation 
techniques can be used to obtain the solution of the spherical (or axisymmetric) problem (e.g. Lynden-Bell 1962; 
Hunter 1975; Kalnajs 1976; Dejonghe 1986; Qian and Hunter 1995) but there is the same difficulty of requiring 
not only the validity of these transformations of the density but also the complex analyticity of a density-related 
integral kernel to complex arguments. The contour integral method of Hunter and Qian (1993) can be also used 
to find anisotropic DFs for spherical systems (Qian and Hunter 1995) but it is valid for densities that are analytic, 
and whose singularities satisfy some conditions. 

It is worth mentioning that Kuzmin and Veltmann (1967a, 1973) and Veltmann (1961, 1965, 1979, 1981) 
developed some more general classes when DF is a product of an unknown function over one argument on a 
known function of another. 

The fundamental integral equations of the problem are given in Sect.|2] and some new formulae for finding 
other anisotropic DFs for stellar systems with known spherically symmetric densities are presented in Sect. [3j 
These depend only on the two integrals of the energy and the magnitude of the angular momentum. This work 
constitutes the core of this paper. These formulae in fact come from an combination of the ideas of Eddington 
and Fricke (see above). Of course, they can be also regarded as simply an extension of Eddington's formula. A 
type of anisotropic DF which is a sum of products of functions only of the energy and powers of the magnitude 
of the angular momentum is derived in Sect. 13.11 and another, which is a sum of products of functions only of 
a special variable and powers of the magnitude of the angular momentum, in Sect. 13.21 More general formulae 
are given in the last part of Sect. [3] Various formulae of the velocity dispersions for such models of these DFs 
are also shown in all the three parts of Sect. [3j Several examples are given in Sect. |4] including some of new 
anisotropic DFs. Different anisotropic DFs for the Plummer model are first given in Sect. 14.1 1 Then anisotropic 
DFs of the Henon isochrone model appear in Sect. 14.21 and those of the y-model are described in Sect. 14.31 
Section[5]is a summary and conclusion. 



2 The fundamental integral equations 

For convenience of calculation, as in BT it is usual to introduce the relative potential \ff and the relative energy 
e of a star in a stellar system, defined by Xjf = —<P + <Pq and e = —E + <Po, where <P and E are, respectively, 
the potential and the energy of a star, and <Pq is a constant generally chosen to be such that there are only stars 
of £ > in the system modelled by DFs. It is well known that the relative energy £ and the three components 
of the angular momentum vector L are four isolating integrals for any orbit in a spherical potential. By the 
Jeans theorem, it follows that the DF of a steady-state spherical stellar system can be regarded as a non-negative 
function of these integrals, denoted by / = /(£, L) . If the system is spherically symmetric in all its properties, / is 
independent of the direction of L and depends only on its magnitude L (e.g., Shiveshwarkar 1936; Ogorodnikov 
1965), that is, / can be expressed as a non-negative function of the relative energy £ and the absolute value 
L of the angular momentum vector L, denoted by / = f{e,L). If the stellar system itself provides the relative 
potential iff = Xff{r) , the mass density p = p (r) can be obtained from Poisson's equation, and it is known that its 
distribution function / = /(r, v) satisfies 



where r is a position vector, v is a velocity vector, G is the gravitational constant. Let r be the modulus of the 
position vector r. For a spherically symmetric system, both the relative potential \ff = \ff{r) and the mass density 
p — p{r) can be, respectively, regarded as two functions of r, that is, iff = iff {r) and p — p{r), and the distribution 
function / = /(r,v) can be expressed as a function of the relative energy £ and the momentum magnitude L, 




(1) 
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thus equation ([Tl can be rewritten as 

Once f{£,L) is known, p can be easily calculated by integration and \ff by solving Poisson's equation for the 
spherical system. The inverse problem that is now investigated is how to derive the DF from the density p for any 
spherical system. Different classes of anisotropic DFs will be below shown, which are derived from spherical 
density profiles for galaxies, by combining some functions only of e (or Q = e — L^/ (2r^)) with some functions 
of the form L^"P" where r^, is a scaling radius, n is an integer greater than —2 and j3„ is a constant such that 
nj3„ > -1. 



3 Anisotropic DFs 

In this section various formulae for anisotropic DFs are obtained from spherical density profiles of different 
forms and their radial and transverse velocity dispersions are expressed in a simple form. We will consider only 
infinite-size models, so e = will be an energy of escape from a system. 



3.1 DFs of the form I L?-"h„{e) 

n=0 



The integral on the right side of Q is first expressed in velocity space. Let v = (vr,V(p,ve) be the velocity in 
the spherical coordinates (r, 9,0) of space and vt = + Vg be the transverse velocity. Then L = rvr and, by 
£ = \jf - {v^ + Vj) /2, ^ can be expressed as 



p(r) =471 / 
Jo 



^ f{\l/-£,rVT)vT 



2e 



dVT 



de 



(3) 



for the spherical density p = p{r) since the system has only stars with e > 0, or say, f{£,L) = for e < 0. This 
means that the mass density p(r) can be regarded as a function depending on the radial coordinate r and the 

m 

relative potential iff. Let p(r) be below denoted by p{iff,r). Assume that f{£,L) = L^^Ki^) for e > and 

n=0 

f{£,L) = for e < 0, and that the density p{\if,r) has the form of p{^f,r) = £ pn{\ff)r-" corresponding to the 

«=o 

relative potential iff. Notice that such representation of the density is not unique. Similar models were discussed 
earlier by Bouvier (1962, 1963). By Q, it follows that 



=L r(« + 3/2) Jo 



hn{£){^f -£)"+' l^d£. 



(4) 



n=0 n=0 

Equation (|4]l holds if /?„ (e) (n = 0, 1 , 2, . . . , m) satisfy the conditions 

Pniw) = (2^)'^ V(«+"3/2) i''^"^^^^'^~^^"^''''^^- 



(5) 



By taking the (n+ l)st derivative of both sides of equation Q with respect to Xjf and using Abel's integral 
equation, it is found that equation ^ can be solved by 



h„{£) = 



'd"+^p„{w) dxff , 1 fd"+'p„{xtf) 



/n+l 



M/=0 



(6) 
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for n = 0, 1,2, . . . ,m. Hence it can be easily found from Q that 



i//=0 



(7) 



for £ > 0, which is a class of anisotropic DFs generated from the spherical density of the form p{\ff,r) = 

m 

L Pn{v)f^"- Furthermore, if {dJ pf,{\if)/d\if')\if=o = for 7 = 0, 1, . . . ,m + 1, then, for e > 0, Gil can be ex- 

n=0 

pressed as 

' ' V8n\h2"n\ de-+^ Jo ' 
By (|7]i, the velocity dispersions (Wt^) ™d ^riV^^) can be also found to be of the following forms 



(8) 



..2n 



PiV,r);f^o Jo 



pn{\j/)d\i/ 



and 



1 rV 



(9) 



(10) 



for any DF derived from the spherical density of the form p[^r,r) = ^ p„[^r)r . It can be also known that 

these dispersions ^ and (fTOl l can be obtained directly according to the following velocity dispersion formulae 
(Dejonghe 1986, 1987) 



1 



p{¥,r)Jo 



p{\l/,r)d\j/ 



and 



a^(v/,/-) = a^(VA,r) 



p{¥,r)Jo 



Vd[r^pixi/,r)] 



d^r'. 



(11) 



(12) 



3.2 DFs of the form I L^"g„{Q) 



n=0 



A more general expression for the integral in the right side of Q can also be derived. To do this, first put 
Q — Iff ~ (^r + ^7-)/2 ~ '■^Vj-/(2rQ), where and vj are the same as in Sect. 13.1 1 and rg is a scaling radius. 
Obviously, Q = e — L^/(2r^) and Q ^ e as ^ 0°. Assume that the DF is of the form / = f{Q,L), and that the 
system has only stars with Q > 0, or equivalently, / = for Q < 0. Then it is easy to see that 



p{r)=4n 

Jo 



Vt 



2G-(l+rVr2)v2 



=dvi 



dQ. 



(13) 



This shows that the density p(r) can be also regarded as a function of the radial coordinate r and the relative 

m 

potential 1//. Suppose that f{Q,L) = L^"8n{Q) for Q > 0, and that the density p{r) can be expressed as 

(1=0 

m 

PiV, r)= p„(v^)r2"/(l + r^/rl)"+K Then similarly by OS it follows that 



|?,(l + '-V^«')"+' tor{n + 3/2){l+ryrlr+^ Jo 



n=0 

PniVV^" 



(27r)3/22«n!r 



2« 



(14) 
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Equation ( fT4] i holds if (2) (« = 0, 1 , 2, . . . , m) satisfy 



1 (n + 3/2) Jo 



(15) 



By taking the (« + l)st derivative of both sides of equation ( 115b with respect to iff, and using Abel's integral 
equation, the solution of equation ( 1151) is given by 



8n{Q) 



d^f"+'^ VQK dxif"+^ 



i//=0 



(16) 



for n = 0,l,2,...,m. Then it can be easily shown from ( 116b that 



m jln 



dY"+^ VQ\ dY 



d"+'Pn{W) 



(17) 



for Q > 0, which is another class of anisotropic DFs generated from the spherical density of the form of p (i/A, r) = 

m 

I PniY)r^"/{i+r^/rl)"^^. Furthermore, if {d-i pniY)/dY')w=o = for j = 0, 1,2, .. . ,m+ 1, then for Q>0, 

n=Q 

(117b can be rewritten as 



f(n n = -L- V f^ PniWw 

^ V^K^ h^"n\ dQ"+^ ]q VQ^' 



(18) 



Of course, ( 117b and dlSb coincide with (|7]l and (HI, respectively. In other words, (|7]l and ([Hi are, respectively, 
limits of (fTTl l and ( fTSl) when r„ ^ oo. 

Similar to those in Sect. 13. ll the velocity dispersions a^{\ff,r) and aj{\lf,r) can be also found to be of the 
following forms 



p{V,r) ^Q{l+r^/ri)"+^ Jo 



and 



2 f2n 



(19) 



(20) 



for any DF derived from the spherical density of the form p(i//',r) = XI Pniw)^'^" / ^ /''a) ■ 

It is worth mentioning that if the sum contains only one term then such models have been studied by Cudde- 
ford (1991). Therefore the above DFs can be indeed regarded as a simple generalization of the models given by 
Cuddeford. However, (117b and (118b are new formulae for our finding anisotropic DFs. 



3.3 Miscellaneous DFs 

One can also obtain more general formulae than (|7]i and ( 117b . It can be shown that 

m 

f{e,L) = l^B„L^"P- 

for £ > is an anisotropic DF for the spherical density of the form p{\ff,r)= Pn{^f)r^"^" , with nji„ > — 1. 

n=0 

Here B„ = [(27r)-'/^2"^"r(nj3„ + l)r(l - a,j)]"\ a„ = nj3„ - a„ + 3/2 and a„ is a non-negative integer such that 
< a„ < 1 for « = 0, 1 , • • • ,m. 



^ d''"+^ p„{\l/) d\i/ 



dY 



■a„ + \ 



1 



d""p„iv) 
dxi/"" 



(21) 



l//=0 
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Similarly, it can be readily shown that a further anisotropic DF is given by 



f{Q,L)=l^B„L^"P" 



n=0 



(22) 



for Q > 0, corresponding to the spherical density p{\lf,r) = £ p„(i//)r^"'^"/(l + / rl)"^"^^ , where Q is given 

in Sect. l3.2l and B„,an,a„ and /3„ are the same as in (1211) . 

Put Q = max(Q, 0) . Then, it can be furthermore shown that the DFs of the form 



n=0 



cE fl'"" + 'p„(V^) d\ii I ( d" 



rf°l>.p„(y) \ 

rf''2B + 'p„(y) rfy, l / rf"2„p„(,j/) \ 

correspond to a spherically symmetric density of the form 

m m 



«=0 



V/=0 



(23) 



(24) 



«=o 



n=0 



with n^in > -1, where S,„ = [(27r)^/^2"^'"r(n/3,„ + l)r(l - a,„)]^', a/„ = n/3,„ - a,,, + 3/2 and a,„ is a non- 
negative integer such that < a,„ < 1 for / = 1,2 and « = 0, 1 , • • • , m. 

Finally, the velocity dispersions a,^(i//,r) and aj{\j/,r) can be also obtained as 



Jo 



pn{w')dw' 



(l + r2/r2)«fo,+i Jo 



Pn{w')dw' 



(25) 



and 



(«/3i„ + l)/-2«/5'" / p„(v/)^v/ + i; 

Jo (1 



(nj32„ + l)r2«ft„ /-v 

' Pn{V)dV 



(1 +r2/r2)"fc„+i ./o 



(26) 



for any DF derived from the spherical density given by (|24| |. 



4 Application to spherical densities 

Some anisotropic DFs for spherically symmetric densities are given by using formulae obtained in Sect. [3] Dif- 
ferent anisotropic Plummer models are given in Sect. 14.11 together with the ratios of their radial and transverse 
velocity dispersions Then anisotropic Henon isochrone models appear in Sect. 14.21 and anisotropic y-models in 
Sect. 14.31 For convenience of calculation, dimensionless quantities are used for all the spherical models consid- 
ered, and it is assumed in this section that Q = e — I? /2 (or say, take = 1 in the previous section). 



4. 1 Anisotropic Plummer models 

The well-known Plummer model is given by the dimensionless potential-density pair: 



V/(r) = l/\/l+r2, (27) 



p(r) = [3/(4;r)]/(l+r2) 



(28) 
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First, DFs of the form given in Sect. 13. H are considered for the Plummer model. Eq. ( 128]) can be expressed as 
p{\\f,r) — [3/(471)])/^+^™ (1 where m is a positive integer. Then, by the anisotropic Plummer model 

is given by 



3r(2m + 6)e 



Im+l/l 



2(2;r)V2 



If 



(2m + 9/2-n)(m-n)!(«!)2 \2e) 



(29) 



for e > 0. It can be easily proved that equation ( |29l l is in agreement with that given by Dejonghe (1986). Figure 
[U shows that the contours given by ( 1291 ) resemble those of the even DFs for some axisymmetric models and that 
the larger the parameter m, the more anisotropic the model. Furthermore, by ^ and dlOb . the velocity dispersion 
ratio a^{xif,r) /aj{\lf, r) for the anisotropic model given by i29i can be obtained as follows: <7,^( V, I'^ji^^^) ~ 
{l+r^)/[\ + {m+\y]. This is a very good analytical proof of the anisotropic property shown by Figure [T] 



Eq. ( 1281) can also be rewritten as p{^i,r 
the anisotropic DF of the Plummer model as 



25/2 



V/^{10(c 



[3/(47r)]i/^(l+cr2)/[v/2^c(l 
l)V + 3c[9(l-c)v/2 + 7c]} 



Y^)] for c> 0, giving by (HI), 



[c+(l-c)v/2]3^/e-nj7 



-dxj/ 



9cL^ i/{10(l -c)V^ + 39c(l -c)2v/*+7c2[5c + 8(l -c)vA^]} 



dxif 



(30) 



for e > 0. Eq. (1301 ) can be calculated analytically in terms of generalized hypergeometric functions. However, it is 
not generally easy to calculate numerically these hypergeometric functions. To evaluate the DFs, it is necessary 
to estimate the integrals in ( l30b . Figure |2] shows the contours of the two different anisotropic DFs. Notice that 
c is an anisotropic parameter ranging from to 1. The larger the parameter c, the more anisotropic the model. 
When c = 0, it degenates to be an isotropic model; when c = 1 , it becomes the same anisotropic model as given 
by ( 129b with the parameter in= I. The velocity dispersion ratio G^{\l/,r)/aj{xi/,r) for the anisotropic model 
defined by ( |30] | can be further given to be of the simple form a^{\j/,r)/aj{\i/,r) = (1 +cr^)/(l +2cr^). 
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Next, the DFs of the form given in Sect. l3.2l are investigated for the Plummer model. The density defined by 
can be expressed as 

(31) 



p(V/,r) = [3v/V(4?r)][c/(l + (1 -c)(l - xifY'r^il+r'f] 
for < c < 1 . By dlSl l. the anisotropic DF for equation ( |3T]) can be obtained as 



f{Q,L) 



8cG 



3/2 



3{\-c)l} 



2 1 +61/^2 + </ d\if 



(32) 



for Q>0. The integral in eq. (132b can be calculated analytically in terms of hypergeometric functions, and for 
any c € [0,1) and L > 0, the anisotropic DF given by ( |32] | increases like ( 1 — Q) as Q ^ 1 • When c = 1 , (|32] | 
is the anisotropic DF of Ossipkov-Merritt type (OM). It was found earlier by Kuzmin and Veltmann (1967a). 
Figure |3] displays the contours of the two different DFs obtained by estimating directly the real integral in ( 1321 ). 
Obviously, all the contours plotted here lie in the area of Q > in the physical domain since the DFs defined by 
(132b are zero out of this area. But, the smaller the anisotropic parameter c, the more anisotropic the model. By 
([T9] l and ( |20] |. the velocity dispersion ratio a} {^if ,r) / aj{'^ ,r) for the anisotropic model defined by (l32l) can be 
also given in the following form 



oKw.r) _ (l+r^){2ci/(l+r^)-(l-c)r^[V/^ + ln(l-V/^)]} 
a2(v/,r) ^ 2{c^i^{\+r^)-{\-cy[\if^+\n{\-\if'^)]} ' 

Inserting (|27] | into ( 1331) gives 

oKr) _ oHw.r) _ (l+'-2){2c-(l-c)r2[l+ln(^)]} 



(33) 



(r) a^{w,r) 



2{c-(l-c)r2[l+ln 



(34) 
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By using the anisotropy parameter defined by Binney (1980), it can be also found from ( 1341) with c = 1 that the 
galaxy is isotropic at the centre and becomes increasingly radially anisotropic with radius and that the velocity 
distribution is arbitrarily close to one made entirely of radial orbits at sufficiently large r (Cuddeford 1991). 



4.2 Anisotropic Henon isochrone models 

Henon's (1959a, b) isochrone model has the dimensionless potential-density pair 

¥{r) = ^ 



p(r) 



1 3{l + Vl+r^) + 2r^ 
471 (l+7T+7^j3(7T+^ ' 



(35) 



(36) 



Some anisotropic DFs for such model were found by Kuzmin and Veltmann (1967b, 1973). Now let us show the 
other anisotropic DFs for the isochrone model. Eq. ( 136b can be expressed as 



471(1 -v/-)3 ' 
thus, by ([8]l, giving the anisotropic DF as 



fi£,L) = 



xj/^ (20- 15x1/ + 6x1/^ -\l/^)d\i/ 



(37) 



(38) 



which can be expressed in terms of elementary functions (See Appendix A). Figure [4^ a) shows the contours of 
the DF given by i 




Fig. 4 The contours of the DFs for Henon's model, (a), (b) and (c) are for l l38t . l l40t and <42t . respectively. In (a), successive contour 
levels differ by factors of 0.2. In (b) and (c), successive contour levels differ by factors of 0.5. 
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Eq. ( 136b can be expressed as 



An 



3v^^ 1 2v/*'r 
1 - V/ 1 + r2 ^ (1 - V/)3 



(39) 



and then it follows from ( 123]) that the anisotropic DF can be obtained as 



f{£,Q,L) 



Q \jf{3-3\if+\if^)d\if 



[l-W)\Q-¥ 



V/^(20- 15X1/ + 6x1/^ -\i/^)d\i/ 



(40) 



where Q = max(Q,0). The integral in ( l40l) can be calculated analytically in terms of elementary functions (See 
Appendix A). Figure |4jb) displays the contours of the DF given by ( 1401) . 
The density ( 136b can be also expressed as 



An 



Ixif^r 



■6^2 



(41) 



and, corresponding to it, by ( 123b . the anisotropic DF is 



f{£,Q,L) 



25/2^:3 



^ \l/^{5-A\i/+\i/^)dii/ 



Jo 



AQ 



3/2 







(42) 



where Q = max(Q,0). The two integrals at the right side of equation (|42] | can be also expressed in terms of 
elementary functions (See Appendix A). FigurelUc) illustrates the contours of the DF given by ( 142b . 



4.3 Anisotropic y-models 

The dimensionless potential-density pair of the y-model (Kuzmin, Veltmann, Tenjes 1986; Dehnen 1993; Saha 
1993; Tremaine et al. 1994) is 

'^^'^^"\ln[(r+l)/r], y=2 ^^^^ 

P(^) = y^ rr(r+l)4-r - ^44) 

In the same way as for the Plummer model, by (l8]l, it can be shown that the anisotropic DFs of the form of 
Sect. l3.1l can be obtained if the density defined by (1441) is expressed as 



3-y(l-y)^+^»'(l+r^)'" 



PiV,r) = ^ \J! , (45) 



where m is a positive integer and y is written as 



^^^,^(2-r)v.]V(-r),,^2, (46) 



In fact, it is known from (1431) that y = r/{r+l). These anisotropic DFs cannot be expressed in terms of elemen- 
tary functions or generalized hypergeometric functions although they can be calculated numerically. 
Another expression of the density for the y-model is 

3-rfl-v)'^+^"'fl + r)2'« 
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O-O O-S 0.4 o.e o.s -1 .o 



£ 

Fig. 5 The contours of the anisotropic Hernquist DFs l |49t . Successive contour levels differ by factors of 0.1. 



where y is given in ( l46l) . It can be shown from ( |2T]) that the anisotropic DFs corresponding to ( |47]) can be 
expressed generally in terms of elementary functions or generalized hypergeometric functions. In particular, 
when 7= 1, it is a model obtained by Kuzmin and Veltmann (1973) and Hernquist (1990) and its density can be 
expressed as 

1 f4-\-2.tn 

For example, by use of ( |2T]) . a DF of the mass density given by ( |48]) with m = 1 is 

for e > 0, where h{\if) = i/^(15 - 241//+ 10v/'^)/(l - Xff)^. The two integrals in eq. (|49] | can be expressed in 
terms of elementary functions (See (157b and ( 158b in Appendix A). The contours of the DF given by ( 149 b are 
plotted in Figure |5] 

By dlSb . the anisotropic DFs of the form in Sect. l3.2l can also be obtained in terms of elementary functions or 
generalized hypergeometric functions in the same way as the DFs of Ossipkov-Merritt type are given by Dehnen 
(1993). The detail derivation of these expressions is omitted here. 

By the way, Baes and Dejonghe (2005) recently used the spherical y-models to investigate the dynamical 
structure of isotropic spherical galaxies with a central black hole. Such work is still very significant. 



■ h{\ff)d\tf 



-2V2Lh{e) + - 



3L^ /-^ V^^(20-45vA + 36v/^ + 10v/^)iiV^ 



(49) 



5 Conclusions 

Although real galaxies are hardly spherically symmetric, it is a very necessary and significant step to study the 
self-consistent anisotropic DFs for stellar systems with the known spherically symmetric density. This is not 
only because some important properties (e.g. surface density and cumulative mass) of the spherical model are 
similar to those of the models that can be generated from the spherically symmetric density by replacing the 
spherical radius r by an axisymmetric or triaxial radius in = \Jx^ + {y/q2Y + (-^/^s)^; but also because both 
spherical and non-spherical models have the same typical behaviours of the dynamical quantities in the limits 
of small and large radii. Such models are elliposoidal. Other classes of flattened models can be constructed by 
using the equipotential method (e.g., Kutusov and Ossipkov 1980; Jiang 2000). 

Anisotropic DFs can be obtained for stellar systems with known spherically symmetric density as a sum of 
products of functions only of the potential and a special function (or power) only of the radial coordinate, i.e. 
these DFs are a sum of products of functions only of a special variable (or the energy) and a power only of the 
magnitude of the angular momentum. This comes from a combination of the ideas of Eddington and Fricke. 
It is an extension of Eddington's classical solution for the isotropic DF of a known spherical density. Like his 
method, it requires that the density be expressed as a special function of the potential, and now also of the radial 
coordinate. Also, part of it is a sum of real integrals of functions only of the potential. Most of these integrals 
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can be calculated analytically in terms of elementary functions or generalized hypergeometric functions (e.g. 
Kuzmin and Veltmann 1967b). 

These formulae concerning anisotropic DFs are obtained by use of the Abel integral equation, and they are 
suitable for all such densities that allow the real integrals of the potential to be valid. The anisotropic DFs of 
Ossipkov-Merritt type are their examples in special cases. After the application to Plummer's spherical model, 
it is found that the anisotropic DFs given by ( |29] ) are in fact the same as those given by Dejonghe (1986). Three 
different anisotropic DFs are given for Henon's (1959a, b) isochrone model and they can be expressed in terms 
of elementary functions. Many anisotropic DFs of y-models, expressed in elementary functions or generalized 
hypergeometric functions, can also be obtained. One can further give formulae of the velocity dispersions for 
these anisotropic DFs. 

These expressions can be extended further to axisymmetric systems. It is straightforward to find their ana- 
logues for axisymmetric systems and they can be used to obtain the even DF of Binney's (BT) logarithmic 
potential although Evans (1993) derived it using Lynden-Bell's (1962) method. For the well-known Lynden-Bell 
(1962) model, these analogues degenerate into the method of Fricke (1952). 
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Appendix A: Explicit expressions of integrals for the anisotropic DFs of both Henon's model and Hern- 
quist's model 

Explicit expressions of integrals for the DFs of the Henon model are given below. Using the identity 
\lf^{lO-5\jf+\if^) 6 15 10 



the expression 

I-" \if^ {10 -5\if+\i/^)d\if _ (19- 106e + 240e2-48e'')ye 



1, (50) 



Jo ~(^-^/)5^/e^ " 32(1 -e)^ 

, 15(3-126 + 16e2)arcsin(v^) ^^^^ 



32(l-e)'5/2 
can be derived. Similarly, the identity 

xj/^ (20 -15x1/ + 6x1/^-11/^) 10 24 15 



(l-V/)" (1-V^)^ (l-xif)^ (l-V^r ' 

gives 

'■^ xff^ilO- I5xif + 6xif^-xif^)dxff (53-330e + 880e2-352e''+64e4)v^ 



/ 

Jo 



[l~Y)^Ve^ 64(1 -e)5 



;i-V/)5 (l-v/)5 (l-v/)4 (1-v/) 



(52) 



, 15(5-24£+40£2)arcsin(V?) 

"I TTT-, viT79 ^v^- '^-'^> 

64(1 -e)"/2 



Note that ^(^-'^+f) = 1 - 1. Then 

fQ xif{3-3xif+xif^)dxif ^ {5-2Q)^ 3arcsin(v^) ^ 

Jo {l-xif)^VQ^ 4(1-2)2 + 4(1-G)V2 ^ ^ 

Using the identity 

i/(5-4v/+V^2) 2 6 5 



3 -2 + (1-1//), (55) 
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the result 



/» 



i/(5-4v/+V^2)iiV^ 



(87 - 350£ + 464e^ - 96e^)y/e 
96(1 -e)4 

5(7-24e + 24e2)arcsin(Ve) 2Ve(3 + 2e) 
^ 32(l-e)9/2 3 



(56) 



is obtained. 

Similarly, integrals for the anisotropic DFs of the Hemquist model can be also expressed as 



^0 




(5-2e)ye 3arcsin(ye) 
4(1 -e)2 + 4(l-e)V2 



(57) 



and 




Jo (l-i//)4Ve^ 



(58) 
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